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Appendix A 

<?xml version = "1.0" encoding = ,, UTF-8" ?> 

- <model name = FieldKorosNoyesModel > 

- <notes> 

5 <hl>Field-Koros-Noyes Model of BZ Reaction </hl> 

- <table border-"©" cellspacing = "0" cellpadding = "2"> 

- <thead> 

- <tr> 

<th align = "left" valign = "middle" bgcolor="#eeeeee">Citation</th> 

10 </tr> 

</thead> 

- <tbody> 

- <tr> 

- <td> 

15 R.J. Field and R.M.Noyes,J.Chem.Phys.60,1877 (1974) ; 

R.J. Field, E.Koros,R.M.Noyes,JACS 94,8649 (1972);R.J. Field, R.M.Noyes,Nature 
237,390 (1972) This implementation is taken manufactured by J.D. Murray, 
"Mathematical Biology" (1989) page 181. 

<a href="" /> 
20 </td> 
</tr> 
</tbody> 
</table> 

- <table border="0" cellspacing = "0" cellpadding-'^'^ 
25 <thead> 

- <tr> 

<th align = "left" valign="middle" bgcolor="#eeeeee">Description</th> 

</tr> 
</thead> 
30 <tbody> 

- <tr> 

<td>Field Noyes Version of Belousov- Zhabotinsky Reaction. Br03 is held 

constant; HOBr is typically ignored, and can be replaced by an empty- set. The 
stoichiometry f is typically taken as 1/ 2 or 1 (denominator 1 or 2 in SBML) 

35 .</td> 
</tr> 
</tbody> 
</table> 

- <table border="0" cellspacing = "0" cellpadding = "2"> 
40 <thead> 

- <tr> 

<th align= M left" valign = "middle" bgcolor="#eeeeee">Rate constant </th> 
<th align="left" valign = "middle" bgcolor="#eeeeee">Reaction</th> 

</tr> 
45 </thead> 

- <tbody> 

- <tr> 

<td>kl = 1.3</td> 

<td>Br + Br03 -> HBr02 + HOBr</td> 

50 </tr> 
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- <tr> 

<td>k2 = 2000000</td> 

<td>Br + HBr02 -> HOBr^2</td> 

</tr> 
5 <tr> 
<td>k3 = 34</td> 

<td>Br03 + HBr02 -> Ce^2 + HBr02 A 2</td> 

</tr> 

- <tr> 

10 <td>k4 = 3000</td> 

<td>HBr02^2 -> Br03 + HOBr</td> 

</tr> 

- <tr> 

<td>k5 = 0.02</td> 
15 <td>Ce -> BrM</td> 

</tr> 

</tbody> 

</table> 

- <table border="0" cellspacing="0" cellpadding = "2"> 

29 <thead> 

- <tr> 

<th align="left" valign="middle" bgcolor="#eeeeee">Variable</th> 
<th align="left" valign="middle" bgcolor="#eeeeee">IC </th> 
<th align = "left" valign="middle" bgcolor="#eeeeee">ODE</th> 

25 </tr> 

</thead> 

- <tbody> 

- <tr> 
<td>Br</td> 

30 <td>0.003</td> 

<td>Br'[t] == -(kl*Br[t]*Br03[t]) + f*k5*Ce[t] - k2*Br[t]*HBr02[t]</td> 
</tr> 

- <tr> 
<td>Ce</td> 

35 <td>0.05</td> 

<td>Ce'[t] == -(k5*Ce[t]) + 2*k3*Br03[t]*HBr02[t]</td> 
</tr> 

- <tr> 
<td>HBr02</td> 

40 <td>0.00K/td> 

<td>HBr02'[t] == kl*Br[t]*Br03[t] - k2*Br[t]*HBr02[t] + 
k3*Br03[t]*HBr02[t] - k4*HBr02[t] /v 2</td> 
</tr> 

- <tr> 

45 <td>HOBr</td> 
<td>0</td> 

<td>HOBr'[t] == kl*Br[t]*Br03[t] + 2*k2*Br[t]*HBr02[t] 
+ k4*HBr02[t]^2</td> 

</tr> 
50 </tbody> 
</table> 
</body> 
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</notes> 

- <MstOfCompartments> 
<compartment name="BZ" /> 

</listOfCompartments> 
§ <listOfSpecies> 

<specie name="Br" initialAmount== M 0.003 , ' compartment="BZ" boundaryCondition = "false" 
/> 

<specie name="Br03" initialAmount="0-l" compartment= n BZ n boundaryCondition="true" 

/> 

10 <specie name="Ce" initialAmount="0.05" compartment="BZ" boundaryCondition = "false" 
/> 

<specie name="HBr02" initialAmount="0-001 M compartment="BZ" 

boundaryCondition="false" /> 
<specie name="HOBr" initialAmount= M 0" compartment= M BZ n boundaryCondition="false" 
15 /> 

</listOfSpecies> 

- <listOfReactions> 

- <reaction name= ,, Reactionl" reversible= ,, false ,, > 

- <MstOfReactants> 

20 <specieReference specie="Br" /> 
<specieReference specie= ,, Br03" /> 
< / 1 i stOf Rea eta n ts > 

- <NstOfProducts> 
<specieReference specie="HBr02 M /> 

25 <specieReference specie="HOBr" /> 
</listOfProducts> 

- <kineticLaw formula="Br*Br03*kl"> 

- <listOfParameters> 

<parameter name="kl" value= ,, 1.3" /> 
30 </MstOfParameters> 
</kineticLaw> 
</reaction> 

- < reaction name="Reaction2" reversible=="false M > 

- <MstOfReactants> 

35 <specieReference specie= ,, Br" /> 
<specieReference specie="HBr02" /> 
</MstOfReactants> 

- <NstOfProducts> 

<specieReference specie= ,, HOBr H stoichiometry="2" /> 
40 </NstOfProducts> 

- <kineticLaw formula = "Br*HBr02*k2"> 

- <listOfParameters> 

<parameter name="k2" va!ue= ,, 2000000" /> 
</NstOfParameters> 
45 </kineticLaw> 
</reaction> 

- <reaction name="Reaction3" reversible="false"> 

- <listOfReactants> 
<specieReference specie="Br03" /> 

50 <specieReference specie="HBr02" /> 
</listOfReactants> 

- <MstOfProducts> 
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<specieReference specie="Ce M stoichiometry="2 n /> 
<specieReference specie = "HBr02" stoichiometry= ,, 2" /> 
</listOfProducts> 

- <kineticLaw formula = ,, Br03*HBr02*k3"> 
§ <listOfParameters> 

<parameter name="k3" value="34" /> 
</HstOfParameters> 
</kineticLaw> 
</reaction> 

1Q <reaction name= ,, Reaction4" reversible="false"> 

- <listOfReactants> 
<specieReference specie = "HBr02 n /> 

< / 1 i stO f Re a eta n t s > 

- <MstOfProducts> 

15 <specieReference specie="Br03" /> 
<specieReference specie="HOBr" /> 
</listOfProducts> 

- <kineticLaw formula="HBr02^2*k4"> 

- <MstOfParameters> 

20 <parameter name="k4 M value= M 3000" /> 
</listOfParameters> 
</kineticLaw> 
</reaction> 

- <reaction name="Reaction5" reversibIe= H false"> 
2§ <HstOfReactants> 

<specieReference specie="Ce" /> 
</listOfReactants> 

- <listOfProducts> 

<specieReference specie="Br" stoichiometry="l M denominator= ,, 2" /> 
30 </listOfProducts> 

- <kineticLaw formula = "Ce*k5"> 

- <listOfParameters> 

<parameter name="k5" value="0.02" /> 
<parameter name=T value="0.5" /> 
35 </listOfParameters> 
</kineticLaw> 
</reaction> 
</listOfReactions> 
</model> 
40 </sbml> 
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